! ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
! 
! The Reaction Rates File
! 
! Generated by KPP-2.2.4_gc symbolic chemistry Kinetics PreProcessor
!       (http://www.cs.vt.edu/~asandu/Software/KPP)
! KPP is distributed under GPL, the general public licence
!       (http://www.gnu.org/copyleft/gpl.html)
! (C) 1995-1997, V. Damian & A. Sandu, CGRER, Univ. Iowa
! (C) 1997-2005, A. Sandu, Michigan Tech, Virginia Tech
!     With important contributions from:
!        M. Damian, Villanova University, USA
!        R. Sander, Max-Planck Institute for Chemistry, Mainz, Germany
! 
! File                 : gckpp_Rates.f90
! Time                 : Mon Jan 11 14:22:52 2021
! Working directory    : /n/home08/kbates/GC/Code.12.3.0_new3_arom/KPP/Tropchem
! Equation file        : gckpp.kpp
! Output root filename : gckpp
! 
! ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~



MODULE gckpp_Rates

  USE gckpp_Parameters
  USE gckpp_Global
  IMPLICIT NONE

CONTAINS



! Begin Rate Law Functions from KPP_HOME/util/UserRateLaws

!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
!  User-defined Rate Law functions
!  Note: the default argument type for rate laws, as read from the equations file, is single precision
!        but all the internal calculations are performed in double precision
!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~

!~~~>  Arrhenius
   REAL(kind=dp) FUNCTION ARR( A0,B0,C0 )
      REAL A0,B0,C0      
      ARR =  DBLE(A0) * EXP(-DBLE(B0)/TEMP) * (TEMP/300.0_dp)**DBLE(C0)
   END FUNCTION ARR        

   REAL(kind=dp) FUNCTION ARRPLUS( A0,B0,C0,D0,E0 )
      REAL A0,B0,C0,D0,E0 
      REAL(kind=dp) K0     
      K0 = DBLE(A0) * (DBLE(D0)+(TEMP*DBLE(E0)))
      ARRPLUS =  K0 * EXP(-DBLE(B0)/TEMP) * (TEMP/300.0_dp)**DBLE(C0)
      IF (ARRPLUS < 0.0) ARRPLUS = 0.0
   END FUNCTION ARRPLUS        

   REAL(kind=dp) FUNCTION TUN( A0,B0,C0 )
      REAL A0,B0,C0      
      TUN =  DBLE(A0) * EXP(DBLE(B0)/TEMP) * EXP(DBLE(C0)/TEMP**3)
   END FUNCTION TUN        

   REAL(kind=dp) FUNCTION TUNPLUS( A0,B0,C0,D0,E0 )
      REAL A0,B0,C0,D0,E0
      REAL(kind=dp) K0
      K0 = DBLE(A0) * (DBLE(D0)+(TEMP*DBLE(E0)))
      TUNPLUS = K0 * EXP(DBLE(B0)/TEMP) * EXP(DBLE(C0)/TEMP**3)
      IF (TUNPLUS < 0.0) TUNPLUS = 0.0
   END FUNCTION TUNPLUS        

   REAL(kind=dp) FUNCTION ISO1( A0,B0,C0,D0,E0,F0,G0 )
      REAL A0,B0,C0,D0,E0,F0,G0
      REAL(kind=dp) K0, K1, K2
      K0 = DBLE(D0)*EXP(DBLE(E0)/TEMP)*EXP(1.E8/TEMP**3)
      K1 = DBLE(F0)*EXP(DBLE(G0)/TEMP)
      K2 = DBLE(C0)*K0/(K0+K1)
      ISO1 =  DBLE(A0) * EXP(DBLE(B0)/TEMP) * (1.-K2)
   END FUNCTION ISO1       

   REAL(kind=dp) FUNCTION ISO2( A0,B0,C0,D0,E0,F0,G0 )
      REAL A0,B0,C0,D0,E0,F0,G0
      REAL(kind=dp) K0, K1, K2
      K0 = DBLE(D0)*EXP(DBLE(E0)/TEMP)*EXP(1.E8/TEMP**3)
      K1 = DBLE(F0)*EXP(DBLE(G0)/TEMP)
      K2 = DBLE(C0)*K0/(K0+K1)
      ISO2 =  DBLE(A0) * EXP(DBLE(B0)/TEMP) * K2
   END FUNCTION ISO2       

   REAL(kind=dp) FUNCTION EPO(A1,E1,M1) 
      REAL A1, E1, M1
      REAL(kind=dp) K1      
      K1 = 1.0_dp/(DBLE(M1) * NUMDEN + 1.0_dp)
      EPO = DBLE(A1) * EXP(DBLE(E1)/TEMP) *  K1
   END FUNCTION EPO

   REAL(kind=dp) FUNCTION VILA ( A0,B0,C0,A1,B1,C1 )
      REAL A0,B0,C0,A1,B1,C1
      REAL(kind=dp) K0, K1
      K0 = DBLE(A0) * EXP(DBLE(B0)/TEMP) * (TEMP/300.0_dp)**DBLE(C0)
      K1 = DBLE(A1) * EXP(DBLE(B1)/TEMP) * (TEMP/300.0_dp)**DBLE(C1)
      K0 = K0*NUMDEN
      VILA = K0*K1*0.35_dp/(K0+K1)
   END FUNCTION VILA

   REAL(kind=dp) FUNCTION PAN ( A0,B0,C0,A1,B1,C1,CF )
      REAL A0,B0,C0,A1,B1,C1,CF
      REAL(kind=dp) K0, K1, KR, NC, F     
      K0 = DBLE(A0) * EXP(DBLE(B0)/TEMP) * (TEMP/300.0_dp)**DBLE(C0)
      K1 = DBLE(A1) * EXP(DBLE(B1)/TEMP) * (TEMP/300.0_dp)**DBLE(C1)
      K0 = K0*NUMDEN
      KR = K0/K1
      NC = 0.75_dp-1.27_dp*(LOG10(DBLE(CF)))
      F  = 10.0_dp**(LOG10(DBLE(CF))/(1+(LOG10(KR)/NC)**2))
      PAN = K0*K1*F/(K0+K1)
   END FUNCTION PAN

   REAL(kind=dp) FUNCTION NIT ( A0,B0,C0,n,X0,Y0)
      REAL A0,B0,C0,n,X0,Y0
      REAL(kind=dp) K0, K1, K2, K3, K4     
      K0 = 2.0E-22_dp * EXP(DBLE(n))
      K1 = 4.3E-1_dp*(TEMP/298.0_dp)**(-8)
      K0 = K0*NUMDEN
      K1 = K0/K1
      K2 = (K0/(1.0_dp+K1))*   &
           4.1E-1_dp**(1.0_dp/(1.0_dp+(LOG10(K1))**2))
      K3 = K2/(K2+DBLE(C0))
      K4 = DBLE(A0)*(DBLE(X0)-TEMP*DBLE(Y0))
      NIT = K4 * EXP(DBLE(B0)/TEMP) * K3
      IF (NIT < 0.0) NIT = 0.0
   END FUNCTION NIT

   REAL(kind=dp) FUNCTION ALK ( A0,B0,C0,n,X0,Y0)
      REAL A0,B0,C0,n,X0,Y0
      REAL(kind=dp) K0, K1, K2, K3, K4  
      K0 = 2.0E-22_dp * EXP(DBLE(n))
      K1 = 4.3E-1_dp*(TEMP/298.0_dp)**(-8)
      K0 = K0*NUMDEN
      K1 = K0/K1
      K2 = (K0/(1.0_dp+K1))*   &
           4.1E-1_dp**(1.0_dp/(1.0_dp+(LOG10(K1))**2))
      K3 = DBLE(C0)/(K2+DBLE(C0))
      K4 = DBLE(A0)*(DBLE(X0)-TEMP*DBLE(Y0))
      ALK = K4 * EXP(DBLE(B0)/TEMP) * K3
      IF (ALK < 0.0) ALK = 0.0
   END FUNCTION ALK

!~~~> Simplified Arrhenius, with two arguments
!~~~> Note: The argument B0 has a changed sign when compared to ARR
   REAL(kind=dp) FUNCTION ARR2( A0,B0 )
      REAL A0,B0
      ARR2 =  DBLE(A0) * EXP( DBLE(B0)/TEMP )              
   END FUNCTION ARR2

   REAL(kind=dp) FUNCTION EP2(A0,C0,A2,C2,A3,C3)
      REAL A0,C0,A2,C2,A3,C3
      REAL(kind=dp) K0,K2,K3
      K0 = DBLE(A0) * EXP(-DBLE(C0)/TEMP)
      K2 = DBLE(A2) * EXP(-DBLE(C2)/TEMP)
      K3 = DBLE(A3) * EXP(-DBLE(C3)/TEMP)
      K3 = K3*NUMDEN*1.0E6_dp
      EP2 = K0 + K3/(1.0_dp+K3/K2 )
   END FUNCTION EP2

   REAL(kind=dp) FUNCTION EP3(A1,C1,A2,C2) 
      REAL A1, C1, A2, C2
      REAL(kind=dp) K1, K2
      K1 = DBLE(A1) * EXP(-DBLE(C1)/TEMP)
      K2 = DBLE(A2) * EXP(-DBLE(C2)/TEMP)
      EP3 = K1 + K2*(1.0E6_dp*NUMDEN)
   END FUNCTION EP3 

   REAL(kind=dp) FUNCTION FALL ( A0,B0,C0,A1,B1,C1,CF)
      REAL A0,B0,C0,A1,B1,C1,CF
      REAL(kind=dp) K0, K1
      K0 = DBLE(A0) * EXP(-DBLE(B0)/TEMP)* (TEMP/300.0_dp)**DBLE(C0)
      K1 = DBLE(A1) * EXP(-DBLE(B1)/TEMP)* (TEMP/300.0_dp)**DBLE(C1)
      K0 = K0*NUMDEN*1.0E6_dp
      K1 = K0/K1
      FALL = (K0/(1.0_dp+K1))*   &
           DBLE(CF)**(1.0_dp/(1.0_dp+(LOG10(K1))**2))
   END FUNCTION FALL

  !---------------------------------------------------------------------------

  ELEMENTAL REAL(kind=dp) FUNCTION k_3rd(temp,cair,k0_300K,n,kinf_300K,m,fc)

    INTRINSIC LOG10

    REAL(kind=dp), INTENT(IN) :: temp      ! temperature [K]
    REAL(kind=dp), INTENT(IN) :: cair      ! air concentration [molecules/cm3]
    REAL, INTENT(IN) :: k0_300K   ! low pressure limit at 300 K
    REAL, INTENT(IN) :: n         ! exponent for low pressure limit
    REAL, INTENT(IN) :: kinf_300K ! high pressure limit at 300 K
    REAL, INTENT(IN) :: m         ! exponent for high pressure limit
    REAL, INTENT(IN) :: fc        ! broadening factor (usually fc=0.6)
    REAL(kind=dp) :: zt_help, k0_T, kinf_T, k_ratio

    zt_help = 300._dp/temp
    k0_T    = k0_300K   * zt_help**(n) * cair ! k_0   at current T
    kinf_T  = kinf_300K * zt_help**(m)        ! k_inf at current T
    k_ratio = k0_T/kinf_T
    k_3rd   = k0_T/(1._dp+k_ratio)*fc**(1._dp/(1._dp+LOG10(k_ratio)**2))

  END FUNCTION k_3rd

  !---------------------------------------------------------------------------

  ELEMENTAL REAL(kind=dp) FUNCTION k_arr (k_298,tdep,temp)
    ! Arrhenius function

    REAL,     INTENT(IN) :: k_298 ! k at T = 298.15K
    REAL,     INTENT(IN) :: tdep  ! temperature dependence
    REAL(kind=dp), INTENT(IN) :: temp  ! temperature

    INTRINSIC EXP

    k_arr = k_298 * EXP(tdep*(1._dp/temp-3.3540E-3_dp)) ! 1/298.15=3.3540e-3

  END FUNCTION k_arr

!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
!  End of User-defined Rate Law functions
!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~

! End Rate Law Functions from KPP_HOME/util/UserRateLaws


! Begin INLINED Rate Law Functions


  REAL(kind=dp) FUNCTION OH_O1D (J, H2O, TEMP, NUMDEN)
  REAL*8 J, H2O, TEMP, NUMDEN
  REAL*8 K1, K2, K3
  REAL*8 N2, O2

  N2=0.79d0
  O2=0.21d0

  K1=1.63d-10*exp(60d0/TEMP)
  K2=2.15d-11*exp(110d0/TEMP)
  K3=3.30d-11*exp(55d0/TEMP)

  OH_O1D=J*K1*H2O/(K1*H2O+K2*N2*NUMDEN+K3*O2*NUMDEN)
  END FUNCTION OH_O1D

  REAL(kind=dp) FUNCTION HO2_H2O(H2O,TEMP)
  REAL*8 TEMP, H2O
  HO2_H2O=1+1.4D-21*H2O*EXP(2200/TEMP)
  END FUNCTION HO2_H2O


  REAL(kind=dp) FUNCTION GCARR( A0,B0,C0 )
      REAL A0,B0,C0 
      GCARR =  DBLE(A0) * EXP(DBLE(C0)/TEMP) * (300._dp/TEMP)**DBLE(B0)
  END FUNCTION GCARR    

  REAL(kind=dp) FUNCTION GC_HO2HO2( A0,B0,C0,A1,B1,C1 )
      REAL A0,B0,C0,A1,B1,C1
      REAL(kind=dp) :: R0,R1
      R0 =  DBLE(A0) * EXP(DBLE(C0)/TEMP) * (300._dp/TEMP)**DBLE(B0)
      R1 =  DBLE(A1) * EXP(DBLE(C1)/TEMP) * (300._dp/TEMP)**DBLE(B1)

      GC_HO2HO2 = (R0+R1*NUMDEN)*(1.D0+1.4E-21_dp*H2O* &
                   EXP(2200.E+0_dp/TEMP))
  END FUNCTION GC_HO2HO2
  
  REAL(kind=dp) FUNCTION GC_TBRANCH( A0,B0,C0,A1,B1,C1 )
! Temperature Dependent Branching Ratio
      REAL A0,B0,C0,A1,B1,C1
      REAL(kind=dp) :: R0,R1
      R0 =  DBLE(A0) * EXP(DBLE(C0)/TEMP) * (300._dp/TEMP)**DBLE(B0)
      R1 =  DBLE(A1) * EXP(DBLE(C1)/TEMP) * (300._dp/TEMP)**DBLE(B1)

      GC_TBRANCH = R0/(1.D0+R1)

  END FUNCTION GC_TBRANCH

  REAL(kind=dp) FUNCTION GC_RO2HO2( A0,B0,C0,A1,B1,C1 )
! Carbon Dependence of RO2+HO2
      REAL A0,B0,C0,A1,B1,C1
      REAL(kind=dp) :: R0,R1
      R0 =  DBLE(A0) * EXP(DBLE(C0)/TEMP) * (300._dp/TEMP)**DBLE(B0)
      R1 =  DBLE(A1) * EXP(DBLE(C1)/TEMP) * (300._dp/TEMP)**DBLE(B1)

      GC_RO2HO2 = R0*(1E0-EXP(-0.245E0*R1))

  END FUNCTION GC_RO2HO2

  REAL(kind=dp) FUNCTION GC_DMSOH( A0,B0,C0,A1,B1,C1 )
! DMS+OH+O2
    REAL A0,B0,C0,A1,B1,C1
    REAL(kind=dp) :: R0,R1
    R0 =  DBLE(A0) * EXP(DBLE(C0)/TEMP) * (300._dp/TEMP)**DBLE(B0)
    R1 =  DBLE(A1) * EXP(DBLE(C1)/TEMP) * (300._dp/TEMP)**DBLE(B1)
!    GC_DMSOH = R0/(1e0_dp+R1*0.2095e0_dp)
    GC_DMSOH = (R0*NUMDEN*0.2095e0_dp)/(1e0_dp+R1*0.2095e0_dp)
    
  END FUNCTION GC_DMSOH

  REAL(kind=dp) FUNCTION GC_GLYXNO3( A0,B0,C0 )
! ---  K = K1*([O2]+3.5D18)/(2*[O2]+3.5D18) --- HO2+2*CO branch of GLYX+OH/NO3
    REAL A0,B0,C0
    REAL(kind=dp) R0
    REAL(kind=dp) O2

    O2 = NUMDEN*0.2095e0_dp
    R0 =  DBLE(A0) * EXP(DBLE(C0)/TEMP) * (300._dp/TEMP)**DBLE(B0)
    GC_GLYXNO3 = R0*(O2+3.5E+18_dp)/(2.E+0_dp*O2+3.5E+18_dp)
    
  END FUNCTION GC_GLYXNO3

  REAL(kind=dp) FUNCTION GC_OHHNO3( A0,B0,C0,A1,B1,C1,A2,B2,C2 )
! ---  OH + HNO3:   K = K0 + K3[M] / (1 + K3[M]/K2)  ------
    REAL A0,B0,C0,A1,B1,C1,A2,B2,C2
    REAL(kind=dp) R0,R1,R2
    R0 =  DBLE(A0) * EXP(DBLE(C0)/TEMP) * (300._dp/TEMP)**DBLE(B0)
    R1 =  DBLE(A1) * EXP(DBLE(C1)/TEMP) * (300._dp/TEMP)**DBLE(B1)
    R2 =  NUMDEN*(DBLE(A2) * EXP(DBLE(C2)/TEMP) * (300._dp/TEMP)**DBLE(B2))
    GC_OHHNO3 = R0 + R2/(1.E0_dp + R2/R1)
    
  END FUNCTION GC_OHHNO3

  REAL(kind=dp) FUNCTION GC_GLYCOHA( A0,B0,C0 )
! 
    REAL A0,B0,C0,R0,GLYC_FRAC
    R0 =  DBLE(A0) * EXP(DBLE(C0)/TEMP) * (300._dp/TEMP)**DBLE(B0)
    GLYC_FRAC=1e+0_dp-11.0729e+0_dp*EXP(-(1._dp/73._dp)*TEMP)
    IF (GLYC_FRAC<0e+0_dp) GLYC_FRAC=0e+0_dp
    GC_GLYCOHA=R0*GLYC_FRAC

  END FUNCTION GC_GLYCOHA

  REAL(kind=dp) FUNCTION GC_GLYCOHB( A0,B0,C0 )
! 
    REAL A0,B0,C0
    REAL(kind=dp) :: R0,GLYC_FRAC
    R0 =  DBLE(A0) * EXP(DBLE(C0)/TEMP) * (300._dp/TEMP)**DBLE(B0)
    GLYC_FRAC=1e+0_dp-11.0729e+0_dp*EXP(-(1._dp/73._dp)*TEMP)
    IF (GLYC_FRAC<0e+0_dp) GLYC_FRAC=0e+0_dp
    GC_GLYCOHB=R0*(1e0_dp-GLYC_FRAC)
    
  END FUNCTION GC_GLYCOHB

  REAL(kind=dp) FUNCTION GC_HACOHA( A0,B0,C0 )
! 
    REAL A0,B0,C0
    REAL(kind=dp) :: R0,HAC_FRAC
    R0 =  DBLE(A0) * EXP(DBLE(C0)/TEMP) * (300._dp/TEMP)**DBLE(B0)
    HAC_FRAC=1e+0_dp-23.7e+0_dp*EXP(-(1._dp/60._dp)*TEMP)
    IF (HAC_FRAC<0e+0_dp) HAC_FRAC=0e+0_dp
    GC_HACOHA=R0*HAC_FRAC

  END FUNCTION GC_HACOHA

  REAL(kind=dp) FUNCTION GC_HACOHB( A0,B0,C0 )
! 
    REAL A0,B0,C0
    REAL(kind=dp) :: R0,HAC_FRAC
    R0 =  DBLE(A0) * EXP(DBLE(C0)/TEMP) * (300._dp/TEMP)**DBLE(B0)
    HAC_FRAC=1e+0_dp-23.7e+0_dp*EXP(-(1._dp/60._dp)*TEMP)
    IF (HAC_FRAC<0e+0_dp) HAC_FRAC=0e+0_dp
    GC_HACOHB=R0*(1.E0_dp-HAC_FRAC)
    
  END FUNCTION GC_HACOHB

  REAL(kind=dp) FUNCTION GC_OHCO( A0,B0,C0 )

    REAL A0,B0,C0,R0
    REAL KLO1,KLO2,KHI1,KHI2,XYRAT1,XYRAT2,BLOG1,BLOG2,FEXP1,FEXP2
    REAL KCO1,KCO2,KCO

    R0 =  DBLE(A0) * EXP(DBLE(C0)/TEMP) * (300._dp/TEMP)**DBLE(B0)
    R0 = R0 * (1.E+0_dp + 0.6e+0_dp*9.871E7_dp*PRESS)

    ! new OH+CO rate from JPL2006.
    KLO1=5.9E-33_dp*(300._dp/TEMP)**(1.E+0_dp) 
    KHI1=1.1E-12_dp*(300._dp/TEMP)**(-1.3E0_dp)
    XYRAT1=KLO1*NUMDEN/KHI1
    BLOG1=LOG10(XYRAT1)
    FEXP1=1.E+0_dp/(1.E+0_dp+BLOG1*BLOG1)
    KCO1=KLO1*NUMDEN*0.6**FEXP1/(1.e+0_dp+XYRAT1)
    KLO2=1.5E-13_dp*(300._dp/TEMP)**(0.E+0_dp)
    KHI2=2.1e+09_dp *(300._dp/TEMP)**(-6.1E+0_dp)
    XYRAT2=KLO2*NUMDEN/KHI2
    BLOG2=LOG10(XYRAT2)
    FEXP2=1.E+0_dp/(1.E+0_dp+BLOG2*BLOG2)
    KCO2=KLO2*0.6**FEXP2/(1.e+0_dp+XYRAT2)
    KCO=KCO1+KCO2
    GC_OHCO=KCO
    
  END FUNCTION GC_OHCO

  REAL(kind=dp) FUNCTION GC_RO2NO( B,A0,B0,C0,A1,B1,C1 )
    ! ---  K = K1*(1-FYRNO3(K2,M,T))  ---  abstraction branch of RO2+NO
    CHARACTER(*) B !Branch Toggle
    REAL A0,B0,C0,A1,B1,C1
    REAL(kind=dp) :: R0,R1
    REAL(kind=dp) :: YYYN, XXYN,  AAA,  RARB, ZZYN
    REAL(kind=dp) :: XF, ALPHA, Y300, BETA, XMINF, XM0
    REAL(kind=dp) :: FYRNO3
    Y300 = 0.826
    ALPHA= 1.94e-22
    BETA = 0.97
    XM0  = 0.
    XMINF= 8.1
    XF   = 0.411
    
    R0 =  DBLE(A0) * EXP(DBLE(C0)/TEMP) * (300._dp/TEMP)**DBLE(B0)
    R1 =  DBLE(A1) * EXP(DBLE(C1)/TEMP) * (300._dp/TEMP)**DBLE(B1)
    
    ! Initialize static variables
    
    XXYN   = ALPHA*EXP(BETA*R1)*NUMDEN*((300./TEMP)**XM0)
    YYYN   = Y300*((300./TEMP)**XMINF)
    AAA    = LOG10(XXYN/YYYN)
    ZZYN   = 1./(1.+ AAA*AAA )
    RARB   = (XXYN/(1.+ (XXYN/YYYN)))*(XF**ZZYN)
    FYRNO3 = RARB/(1. + RARB)
    IF (trim(B) .eq. 'A') THEN
       GC_RO2NO     = R0 * FYRNO3
    ELSEIF (trim(B) .eq. 'B') THEN
       GC_RO2NO     = R0 * (1.E+0_dp - FYRNO3)
    ENDIF

  END FUNCTION GC_RO2NO

  REAL(kind=dp) FUNCTION GCJPL3( k0_300, n, ki_300, m)
!  Functions given in JPL Booklet
      REAL k0_300, n, ki_300,m
      REAL k0, ki

      k0=k0_300*((TEMP/300.d0)**(-n))
      ki=ki_300*((TEMP/300.d0)**(-m))

!      GCJPL3=(k0*NUMDEN)/(1+k0*NUMDEN/ki)*0.6** &
!	((1+((LOG10(k0*NUMDEN/ki))**2d0)**-1.0d0))
      GCJPL3=(k0/(1.d0+k0/(ki/NUMDEN)))*0.6** &
        ((1+((log10(k0/(ki/NUMDEN)))**2d0)**1.0e0))
      GCJPL3=GCJPL3*NUMDEN
  END FUNCTION GCJPL3
  
  REAL(kind=dp) FUNCTION GCJPLEQ( A0,B0,C0,A1,B1,C1,A2,B2,C2,FV,FCT1,FCT2 )
 ! Function calculates the rate constant of the forward reaction
 ! calculates the equilibrium constant 
 ! Find the backwards reaction by K=kforward/kbackwards
       REAL A0,B0,C0,A1,B1,C1
       REAL(kind=dp) :: R0,R1
       REAL, OPTIONAL :: A2,B2,C2,FV,FCT1,FCT2 !If a P-dependent rxn

       ! Calculate Backwards reaction
       R0 = GCARR( A0,B0,C0 )

       ! Calculate forwards reaction
       IF (present(A2)) THEN ! P-dependent
          IF (present(B2) .and. present(C2) .and. present(FV) &
                          .and. present(FCT1) .and. present(FCT2)) THEN
	     R1 = GCJPLPR( A1,B1,C1,A2,B2,C2,FV,FCT1,FCT2)
          ELSE ! Missing params!
	     write(*,'(a)') 'GCJPLEQ: Missing parameters for P-dependent reaction.'
	     write(*,'(a)') 'GCJPLEQ: Returning zero'
	     GCJPLEQ = 0.E0
	     RETURN
	  ENDIF
       ELSE
          R1 = gcarr( A1,B1,C1 ) !Std. Arrhenius eqn.
       ENDIF

       GCJPLEQ=R1/R0
  END FUNCTION GCJPLEQ
 
  REAL(kind=dp) FUNCTION GCJPLPR(A0,B0,C0,A1,B1,C1,FV,FCT1,FCT2)
! * PRESSURE-DEPENDENT EFFECTS 
! * ADD THE THIRD BODY EFFECT FOR PRESSURE DEPENDENCE OF RATE
! * COEFFICIENTS.
! A0 B0, & C0 are the Arrhenius parameters for the lower-limit
! rate. A1, B1 & C1 are the upper-limit parameters.
! FV is the falloff curve paramter, (SEE ATKINSON ET. AL (1992)
! J. PHYS. CHEM. REF. DATA 21, P. 1145). USUALLY = 0.6 
! 
       REAL A0,B0,C0,A1,B1,C1,FV,FCT1,FCT2
       REAL FCT,XYRAT,BLOG,RLOW,RHIGH,FEXP
       
       RLOW  = GCARR( A0,B0,C0 )*NUMDEN
       RHIGH = GCARR( A1,B1,C1 )

       IF     (FCT2.NE.0.) THEN	
             FCT            = EXP(-TEMP / FCT1) + EXP(-FCT2 / TEMP) 
             XYRAT          = RLOW/RHIGH
             BLOG           = LOG10(XYRAT)
             FEXP           = 1.e+0_dp / (1.e+0_dp + BLOG * BLOG)
             GCJPLPR        = RLOW*FCT**FEXP/(1e+0_dp+XYRAT) 
       ELSEIF (FCT1.NE.0.) THEN 
             FCT            = EXP(-TEMP / FCT1)
             XYRAT          = RLOW/RHIGH
             BLOG           = LOG10(XYRAT)
             FEXP           = 1.e+0_dp / (1.e+0_dp + BLOG * BLOG)
             GCJPLPR        = RLOW*FCT**FEXP/(1e+0_dp+XYRAT)  
       ELSE
             XYRAT          = RLOW/RHIGH
             BLOG           = LOG10(XYRAT)
             FEXP           = 1.e+0_dp / (1.e+0_dp + BLOG * BLOG)
             GCJPLPR        = RLOW*FV**FEXP/(1e+0_dp+XYRAT)
       ENDIF
       

  END FUNCTION GCJPLPR

  REAL(kind=dp) FUNCTION GCIUPAC3(ko_300,n,ki_300,m,Fc) 
! Function calcualtes the rate constant of 3 body reaction using IUPAC 
! methology
  REAL ko_300,n,ki_300,m,Fc
  REAL ko, ki, F, NN

  ko=ko_300*((TEMP/300.e0)**n)*NUMDEN
  ki=ki_300*((TEMP/300.e0)**m)

  NN=0.75-1.27*LOG10(Fc)
  F=10.0**(LOG10(Fc)/(1.0e0+(LOG10(ko/ki)/NN)**2.0))

  GCIUPAC3=ko/(1+ko/ki)*F
  END FUNCTION GCIUPAC3


! End INLINED Rate Law Functions

! ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
! 
! Update_SUN - update SUN light using TIME
!   Arguments :
! 
! ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~

  SUBROUTINE Update_SUN()
      !USE gckpp_Parameters
      !USE gckpp_Global

    IMPLICIT NONE

    REAL(kind=dp) :: SunRise, SunSet
    REAL(kind=dp) :: Thour, Tlocal, Ttmp 
    ! PI - Value of pi
    REAL(kind=dp), PARAMETER :: PI = 3.14159265358979d0
    
    SunRise = 4.5_dp 
    SunSet  = 19.5_dp 
    Thour = TIME/3600.0_dp 
    Tlocal = Thour - (INT(Thour)/24)*24

    IF ((Tlocal>=SunRise).AND.(Tlocal<=SunSet)) THEN
       Ttmp = (2.0*Tlocal-SunRise-SunSet)/(SunSet-SunRise)
       IF (Ttmp.GT.0) THEN
          Ttmp =  Ttmp*Ttmp
       ELSE
          Ttmp = -Ttmp*Ttmp
       END IF
       SUN = ( 1.0_dp + COS(PI*Ttmp) )/2.0_dp 
    ELSE
       SUN = 0.0_dp 
    END IF

 END SUBROUTINE Update_SUN

! End of Update_SUN function
! ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~


! ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
! 
! Update_RCONST - function to update rate constants
!   Arguments :
! 
! ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~

SUBROUTINE Update_RCONST ( )




! Begin INLINED RCONST


! End INLINED RCONST

  RCONST(1) = (GCARR(3.00E-12,0.0E+00,-1500.0))
  RCONST(2) = (GCARR(1.70E-12,0.0E+00,-940.0))
  RCONST(3) = (GCARR(1.00E-14,0.0E+00,-490.0))
  RCONST(4) = (GCARR(1.20E-13,0.0E+00,-2450.0))
  RCONST(5) = (GCARR(2.90E-16,0.0E+00,-1000.0))
  RCONST(6) = (GCARR(1.80E-12,0.0E+00,0.0))
  RCONST(7) = (GCJPLPR(6.90E-31,1.0E+00,0.0,2.6E-11,0.0,0.0,0.6,0.0,0.0))
  RCONST(8) = (GCARR(4.80E-11,0.0E+00,250.0))
  RCONST(9) = (GCARR(1.80E-12,0.0E+00,0.0))
  RCONST(10) = (GCARR(3.30E-12,0.0E+00,270.0))
  RCONST(11) = (GC_HO2HO2(3.00E-13,0.0E+00,460.0,2.1E-33,0.0,920.0))
  RCONST(12) = (GCARR(2.80E-12,0.0E+00,-1800.0))
  RCONST(13) = (GC_OHCO(1.50E-13,0.0E+00,0.0))
  RCONST(14) = (GCARR(2.45E-12,0.0E+00,-1775.0))
  RCONST(15) = (GCARR(2.80E-12,0.0E+00,300.0))
  RCONST(16) = (GCARR(4.10E-13,0.0E+00,750.0))
  RCONST(17) = (GC_TBRANCH(9.50E-14,0.0E+00,390.0,2.62E1,0.0,-1130.0))
  RCONST(18) = (GC_TBRANCH(9.50E-14,0.0E+00,390.0,4.E-2,0.0,1130.0))
  RCONST(19) = (GCARR(2.66E-12,0.0E+00,200.0))
  RCONST(20) = (GCARR(1.14E-12,0.0E+00,200.0))
  RCONST(21) = (GCARR(2.66E-12,0.0E+00,200.0))
  RCONST(22) = (GCARR(1.14E-12,0.0E+00,200.0))
  RCONST(23) = (GCARR(5.50E-12,0.0E+00,125.0))
  RCONST(24) = (GCJPLPR(1.80E-30,3.0E+00,0.0,2.8E-11,0.0,0.0,0.6,0.0,0.0))
  RCONST(25) = (GC_OHHNO3(2.41E-14,0.0E+00,460.0,2.69E-17,0.E0,2199.,6.51E-34,0.E0,1335.0))
  RCONST(26) = (GCJPLPR(7.00E-31,2.6E+00,0.0,3.60E-11,0.1,0.0,0.6,0.0,0.0))
  RCONST(27) = (GCARR(1.80E-11,0.0E+00,-390.0))
  RCONST(28) = (GCJPLPR(1.90E-31,3.4E+00,0.0,4.0e-12,0.3,0.0,0.6,0.0,0.0))
  RCONST(29) = (GCJPLPR(9.05E-05,3.4E+00,-10900.0,1.90E15,0.3,-10900.0,0.6,0.0,0.0))
  RCONST(30) = (GCARR(1.30E-12,0.0E+00,380.0))
  RCONST(31) = (GCARR(3.50E-12,0.0E+00,0.0))
  RCONST(32) = (GCARR(1.50E-11,0.0E+00,170.0))
  RCONST(33) = (GCARR(2.20E-11,0.0E+00,0.0))
  RCONST(34) = (GCJPLPR(2.40E-30,3.0E+00,0.0,1.6E-12,-0.1,0.0,0.6,0.0,0.0))
  RCONST(35) = (GCJPLPR(4.14E-04,3.0E+00,-10840.0,2.76E14,-0.1,-10840.0,0.6,0.0,0.0))
  RCONST(36) = (GCARR(4.00E-13,0.0E+00,0.0))
  RCONST(37) = (GCARR(2.90E-12,0.0E+00,-345.0))
  RCONST(38) = (GCARR(4.50E-14,0.0E+00,-1260.0))
  RCONST(39) = (GCARR(5.80E-16,0.0E+00,0.0))
  RCONST(40) = (GCARR(4.63E-12,0.0E+00,350.0))
  RCONST(41) = (GCARR(1.40E-12,0.0E+00,-1900.0))
  RCONST(42) = (GCJPLPR(9.70E-29,5.6E+00,0.0,9.3E-12,1.5E0,0.0,0.6,0.0,0.0))
  RCONST(43) = (GCJPLEQ(9.30E-29,0.0E+00,14000.0,9.7E-29,5.6E0,0.0,9.3E-12,1.5E0,0.,0.6,0.,0.))
  RCONST(44) = (GCARR(8.10E-12,0.0E+00,270.0))
  RCONST(45) = (GCARR(7.66E-12,0.0E+00,-1020.0))
  RCONST(46) = (GCARR(2.60E-12,0.0E+00,365.0))
  RCONST(47) = (GC_TBRANCH(7.60E-12,0.0E+00,-585.0,5.87E0,0.64E0,-816.0))
  RCONST(48) = (GC_TBRANCH(7.60E-12,0.0E+00,-585.0,1.7E-1,-0.64E0,816.0))
  RCONST(49) = (GCARR(2.90E-12,0.0E+00,350.0))
  RCONST(50) = (GCARR(2.70E-12,0.0E+00,350.0))
  RCONST(51) = (GCARR(9.10E-12,0.0E+00,-405.0))
  RCONST(52) = (GC_RO2NO('B',2.70E-12,0.0E+00,350.0,4.5E0,0.0,0.0))
  RCONST(53) = (GC_RO2NO('A',2.70E-12,0.0E+00,350.0,4.5E0,0.0,0.0))
  RCONST(54) = (GCARR(2.70E-12,0.0E+00,350.0))
  RCONST(55) = (GCARR(2.80E-12,0.0E+00,300.0))
  RCONST(56) = (GCARR(2.70E-12,0.0E+00,350.0))
  RCONST(57) = (GCARR(2.70E-12,0.0E+00,350.0))
  RCONST(58) = (GCARR(4.07E+08,0.0E+00,-7694.0))
  RCONST(59) = (GCARR(2.70E-12,0.0E+00,350.0))
  RCONST(60) = (GCARR(2.70E-12,0.0E+00,350.0))
  RCONST(61) = (GCARR(2.90E+07,0.0E+00,-5297.0))
  RCONST(62) = (GCARR(2.70E-12,0.0E+00,350.0))
  RCONST(63) = (GCARR(2.70E-12,0.0E+00,350.0))
  RCONST(64) = (GCARR(2.70E-12,0.0E+00,350.0))
  RCONST(65) = (GCARR(2.30E-12,0.0E+00,0.0))
  RCONST(66) = (GCARR(3.15E-13,0.0E+00,-448.0))
  RCONST(67) = (GCARR(4.00E-12,0.0E+00,0.0))
  RCONST(68) = (GCARR(6.70E-12,0.0E+00,340.0))
  RCONST(69) = (GCJPLPR(9.00E-28,8.9E+00,0.0,7.7E-12,0.2,0.0,0.6,0.0,0.0))
  RCONST(70) = (GCARR(5.20E-13,0.0E+00,980.0))
  RCONST(71) = (GCARR(2.30E-12,0.0E+00,0.0))
  RCONST(72) = (GCARR(2.60E-12,0.0E+00,380.0))
  RCONST(73) = (GCARR(2.06E-13,0.0E+00,1300.0))
  RCONST(74) = (GCARR(2.00E-13,0.0E+00,0.0))
  RCONST(75) = (GCARR(4.15E-15,0.0E+00,-1520.0))
  RCONST(76) = (GCARR(7.48E-12,0.0E+00,410.0))
  RCONST(77) = (GCARR(2.60E-12,0.0E+00,380.0))
  RCONST(78) = (GCARR(2.06E-13,0.0E+00,1300.0))
  RCONST(79) = (GCARR(2.06E-13,0.0E+00,0.0))
  RCONST(80) = (GCARR(1.20E-12,0.0E+00,0.0))
  RCONST(81) = (GCARR(2.70E-12,0.0E+00,350.0))
  RCONST(82) = (GCARR(2.80E-12,0.0E+00,-3280.0))
  RCONST(83) = (GCARR(1.60E-12,0.0E+00,0.0))
  RCONST(84) = (GCARR(3.15E-14,0.0E+00,920.0))
  RCONST(85) = (GCARR(6.00E-12,0.0E+00,410.0))
  RCONST(86) = (GCJPLPR(9.00E-28,8.9E+00,0.0,7.7E-12,0.2,0.0,0.6,0.0,0.0))
  RCONST(87) = (GCJPLEQ(9.00E-29,0.0E+00,14000.0,9.00E-28,8.9E0,0.0,7.7E-12,0.2,0.,0.6,0.,0.))
  RCONST(88) = (GCJPLPR(9.00E-28,8.9E+00,0.0,7.7E-12,0.2,0.0,0.6,0.0,0.0))
  RCONST(89) = (GCJPLPR(9.00E-28,8.9E+00,0.0,7.7E-12,0.2,0.0,0.6,0.0,0.0))
  RCONST(90) = (GCJPLEQ(9.00E-29,0.0E+00,14000.0,9.00E-28,8.9E0,0.0,7.7E-12,0.2,0.,0.6,0.,0.))
  RCONST(91) = (GCJPLEQ(9.00E-29,0.0E+00,14000.0,9.00E-28,8.9E0,0.0,7.7E-12,0.2,0.,0.6,0.,0.))
  RCONST(92) = (GCARR(6.70E-12,0.0E+00,340.0))
  RCONST(93) = (GCARR(6.70E-12,0.0E+00,340.0))
  RCONST(94) = (GCARR(6.50E-15,0.0E+00,0.0))
  RCONST(95) = (1.33E-13+3.82E-11*exp(-2000.0_dp/TEMP))
  RCONST(96) = (GCARR(5.92E-13,0.0E+00,0.0))
  RCONST(97) = (GCARR(5.92E-13,0.0E+00,0.0))
  RCONST(98) = (GCARR(7.40E-13,0.0E+00,700.0))
  RCONST(99) = (GCARR(7.40E-13,0.0E+00,700.0))
  RCONST(100) = (GCARR(8.60E-13,0.0E+00,700.0))
  RCONST(101) = (GC_RO2HO2(2.91E-13,0.0E+00,1300.0,4.0,0.0,0.0))
  RCONST(102) = (GCARR(2.06E-13,0.0E+00,1300.0))
  RCONST(103) = (GC_RO2HO2(2.91E-13,0.0E+00,1300.0,4.0,0.0,0.0))
  RCONST(104) = (GC_RO2HO2(2.91E-13,0.0E+00,1300.0,4.0,0.0,0.0))
  RCONST(105) = (GC_RO2HO2(2.91E-13,0.0E+00,1300.0,4.0,0.0,0.0))
  RCONST(106) = (GC_RO2HO2(2.91E-13,0.0E+00,1300.0,3.0,0.0,0.0))
  RCONST(107) = (GC_RO2HO2(2.91E-13,0.0E+00,1300.0,5.0,0.0,0.0))
  RCONST(108) = (GC_RO2HO2(2.91E-13,0.0E+00,1300.0,3.0,0.0,0.0))
  RCONST(109) = (GCARR(1.30E-12,0.0E+00,-25.0))
  RCONST(110) = (GCARR(3.00E-13,0.0E+00,0.0))
  RCONST(111) = (GCARR(8.00E-16,0.0E+00,0.0))
  RCONST(112) = (GCARR(8.37E-14,0.0E+00,0.0))
  RCONST(113) = (GCARR(8.37E-14,0.0E+00,0.0))
  RCONST(114) = (GCARR(7.50E-13,0.0E+00,500.0))
  RCONST(115) = (GCARR(8.37E-14,0.0E+00,0.0))
  RCONST(116) = (GCARR(8.37E-14,0.0E+00,0.0))
  RCONST(117) = (GCARR(2.30E-12,0.0E+00,0.0))
  RCONST(118) = (GCARR(8.37E-14,0.0E+00,0.0))
  RCONST(119) = (GCARR(8.37E-14,0.0E+00,0.0))
  RCONST(120) = (GCARR(8.37E-14,0.0E+00,0.0))
  RCONST(121) = (GCARR(8.37E-14,0.0E+00,0.0))
  RCONST(122) = (GCARR(8.37E-14,0.0E+00,0.0))
  RCONST(123) = (GCARR(1.30E-12,0.0E+00,0.0))
  RCONST(124) = (GCARR(8.37E-14,0.0E+00,0.0))
  RCONST(125) = (GCARR(3.35E-12,0.0E+00,0.0))
  RCONST(126) = (GCARR(4.60E-12,0.0E+00,70.0))
  RCONST(127) = (GCARR(4.10E-14,0.0E+00,0.0))
  RCONST(128) = (GCARR(2.70E-14,0.0E+00,0.0))
  RCONST(129) = (GCARR(7.40E-13,0.0E+00,700.0))
  RCONST(130) = (GC_RO2HO2(2.91E-13,0.0E+00,1300.0,3.0,0.0,0.0))
  RCONST(131) = (GC_RO2HO2(2.91E-13,0.0E+00,1300.0,3.0,0.0,0.0))
  RCONST(132) = (GCARR(5.20E-13,0.0E+00,980.0))
  RCONST(133) = (GCARR(4.30E-13,0.0E+00,1040.0))
  RCONST(134) = (GCARR(4.30E-13,0.0E+00,1040.0))
  RCONST(135) = (GCJPLPR(4.60E-27,4.0E+00,0.0,2.6E-11,1.3,0.0,0.5,0.0,0.0))
  RCONST(136) = (GCARR(5.50E-15,0.0E+00,-1880.0))
  RCONST(137) = (GCARR(2.90E-11,0.0E+00,0.0))
  RCONST(138) = (GCARR(3.00E-11,0.0E+00,0.0))
  RCONST(139) = (GCARR(8.20E-18,0.0E+00,0.0))
  RCONST(140) = (GCARR(3.48E-12,0.0E+00,0.0))
  RCONST(141) = (GC_GLYCOHA(8.00E-12,0.0E+00,0.0))
  RCONST(142) = (GC_GLYCOHB(8.00E-12,0.0E+00,0.0))
  RCONST(143) = (GCARR(4.59E-13,0.0E+00,-1156.0))
  RCONST(144) = (GCARR(3.10E-12,0.0E+00,340.0))
  RCONST(145) = (GCARR(1.50E-11,0.0E+00,0.0))
  RCONST(146) = (GC_GLYXNO3(1.40E-12,0.0E+00,-1860.0))
  RCONST(147) = (GCARR(3.36E-12,0.0E+00,-1860.0))
  RCONST(148) = (GCARR(3.00E-11,0.0E+00,360.0))
  RCONST(149) = (GCARR(2.70E-12,0.0E+00,580.0))
  RCONST(150) = (GCARR(9.60E-12,0.0E+00,360.0))
  RCONST(151) = (GC_HACOHA(2.15E-12,0.0E+00,305.0))
  RCONST(152) = (GC_HACOHB(2.15E-12,0.0E+00,305.0))
  RCONST(153) = (GCARR(1.68E-12,0.0E+00,500.0))
  RCONST(154) = (GCARR(1.68E-12,0.0E+00,500.0))
  RCONST(155) = (GCARR(1.87E-13,0.0E+00,500.0))
  RCONST(156) = (GCARR(1.87E-13,0.0E+00,500.0))
  RCONST(157) = (GCARR(1.00E-14,0.0E+00,-1970.0))
  RCONST(158) = (GCARR(8.5E-16,0.0E+00,-1520.0))
  RCONST(159) = (GCARR(1.50E-15,0.0E+00,-2110.0))
  RCONST(160) = (GCARR(3.50E-12,0.0E+00,-450.0))
  RCONST(161) = (GCARR(3.40E-15,0.0E+00,0.0))
  RCONST(162) = (GCARR(1.50E-12,0.0E+00,-1815.0))
  RCONST(163) = (GCARR(1.68E-12,0.0E+00,500.0))
  RCONST(164) = (GCARR(1.68E-12,0.0E+00,500.0))
  RCONST(165) = (GCARR(1.87E-13,0.0E+00,500.0))
  RCONST(166) = (GCARR(1.87E-13,0.0E+00,500.0))
  RCONST(167) = (GCARR(1.90E-11,0.0E+00,390.0))
  RCONST(168) = (GCARR(5.18E-12,0.0E+00,200.0))
  RCONST(169) = (GCARR(8.78E-12,0.0E+00,200.0))
  RCONST(170) = (GCARR(5.18E-12,0.0E+00,200.0))
  RCONST(171) = (GCARR(5.18E-12,0.0E+00,200.0))
  RCONST(172) = (GCARR(8.78E-12,0.0E+00,200.0))
  RCONST(173) = (GCARR(8.78E-12,0.0E+00,200.0))
  RCONST(174) = (GCARR(6.13E-13,0.0E+00,200.0))
  RCONST(175) = (GCARR(8.78E-12,0.0E+00,200.0))
  RCONST(176) = (GCARR(6.13E-12,0.0E+00,200.0))
  RCONST(177) = (GCARR(4.14E-12,0.0E+00,200.0))
  RCONST(178) = (GCARR(5.11E-12,0.0E+00,200.0))
  RCONST(179) = (GCARR(1.70E-11,0.0E+00,390.0))
  RCONST(180) = (GCARR(2.97E-11,0.0E+00,390.0))
  RCONST(181) = (GCARR(2.92E-11,0.0E+00,390.0))
  RCONST(182) = (GCARR(4.82E-11,0.0E+00,-400.0))
  RCONST(183) = (GCARR(3.73E-11,0.0E+00,-400.0))
  RCONST(184) = (GCARR(5.79E-11,0.0E+00,-400.0))
  RCONST(185) = (GCARR(3.20E-11,0.0E+00,-400.0))
  RCONST(186) = (GCARR(2.06E-13,0.0E+00,1300.0))
  RCONST(187) = (GCARR(2.70E-12,0.0E+00,350.0))
  RCONST(188) = (GCARR(5.31E-12,0.0E+00,200.0))
  RCONST(189) = (GCARR(8.78E-12,0.0E+00,200.0))
  RCONST(190) = (GCARR(1.84E-12,0.0E+00,200.0))
  RCONST(191) = (GCARR(4.40E-12,0.0E+00,380.0))
  RCONST(192) = (GCARR(6.13E-13,0.0E+00,200.0))
  RCONST(193) = (GCARR(3.60E-12,0.0E+00,380.0))
  RCONST(194) = (GCARR(1.68E-12,0.0E+00,500.0))
  RCONST(195) = (GCARR(1.87E-13,0.0E+00,500.0))
  RCONST(196) = (GCARR(8.37E-14,0.0E+00,0.0))
  RCONST(197) = (GCARR(8.37E-14,0.0E+00,0.0))
  RCONST(198) = (GC_RO2HO2(2.91E-13,0.0E+00,1300.0,4.0,0.0,0.0))
  RCONST(199) = (GC_RO2NO('B',2.70E-12,0.0E+00,350.0,4.0,0.0,0.0))
  RCONST(200) = (GC_RO2NO('A',2.70E-12,0.0E+00,350.0,4.0,0.0,0.0))
  RCONST(201) = (GCARR(6.13E-13,0.0E+00,200.0))
  RCONST(202) = (GCARR(1.40E-18,0.0E+00,0.0))
  RCONST(203) = (GCARR(6.16E-15,0.0E+00,-1814.0))
  RCONST(204) = (GCARR(2.50E-12,0.0E+00,500.0))
  RCONST(205) = (GCARR(1.80E-12,0.0E+00,500.0))
  RCONST(206) = (GCARR(2.00E-13,0.0E+00,500.0))
  RCONST(207) = (GCARR(1.68E-12,0.0E+00,500.0))
  RCONST(208) = (GCARR(1.68E-12,0.0E+00,500.0))
  RCONST(209) = (GCARR(1.68E-12,0.0E+00,500.0))
  RCONST(210) = (GCARR(1.68E-12,0.0E+00,500.0))
  RCONST(211) = (GCARR(1.68E-12,0.0E+00,500.0))
  RCONST(212) = (GCARR(1.68E-12,0.0E+00,500.0))
  RCONST(213) = (GCARR(1.68E-12,0.0E+00,500.0))
  RCONST(214) = (GCARR(1.68E-12,0.0E+00,500.0))
  RCONST(215) = (GCARR(1.68E-12,0.0E+00,500.0))
  RCONST(216) = (GCARR(1.68E-12,0.0E+00,500.0))
  RCONST(217) = (GCARR(1.68E-12,0.0E+00,500.0))
  RCONST(218) = (GCARR(1.68E-12,0.0E+00,500.0))
  RCONST(219) = (GCARR(1.87E-13,0.0E+00,500.0))
  RCONST(220) = (GCARR(1.87E-13,0.0E+00,500.0))
  RCONST(221) = (GCARR(1.87E-13,0.0E+00,500.0))
  RCONST(222) = (GCARR(1.87E-13,0.0E+00,500.0))
  RCONST(223) = (GCARR(1.87E-13,0.0E+00,500.0))
  RCONST(224) = (GCARR(1.87E-13,0.0E+00,500.0))
  RCONST(225) = (GCARR(1.87E-13,0.0E+00,500.0))
  RCONST(226) = (GCARR(1.87E-13,0.0E+00,500.0))
  RCONST(227) = (GCARR(1.87E-13,0.0E+00,500.0))
  RCONST(228) = (GCARR(1.87E-13,0.0E+00,500.0))
  RCONST(229) = (GCARR(1.87E-13,0.0E+00,500.0))
  RCONST(230) = (GCARR(1.87E-13,0.0E+00,500.0))
  RCONST(231) = (GCARR(1.68E-12,0.0E+00,500.0))
  RCONST(232) = (GCARR(1.87E-13,0.0E+00,500.0))
  RCONST(233) = (GCARR(2.50E-12,0.0E+00,500.0))
  RCONST(234) = (GCARR(2.50E-12,0.0E+00,500.0))
  RCONST(235) = (GCARR(8.50E-13,0.0E+00,-2450.0))
  RCONST(236) = (GCJPLPR(1.00E-30,4.8E+00,0.0,7.2E-12,2.1E0,0.0,0.6,0.0,0.0))
  RCONST(237) = (GCJPLPR(1.05E-02,4.8E+00,-11234.0,7.58E16,2.1E0,-11234.0,0.6,0.0,0.0))
  RCONST(238) = (GCARR(3.70E-19,0.0E+00,0.0))
  RCONST(239) = (GCARR(2.90E-17,0.0E+00,0.0))
  RCONST(240) = (GCARR(1.20E-11,0.0E+00,-280.0))
  RCONST(241) = (GC_DMSOH(8.20E-39,0.0E+00,5376.0,1.05E-5,0.0,3644.0))
  RCONST(242) = (GCARR(1.90E-13,0.0E+00,530.0))
  RCONST(243) = (GCJPLPR(3.30E-31,4.3E+00,0.0,1.6E-12,0.0,0.0,0.6,0.0,0.0))
  RCONST(244) = (GCARR(1.60E-11,0.0E+00,-780.0))
  RCONST(245) = (GCARR(4.50E-12,0.0E+00,460.0))
  RCONST(246) = (GCARR(4.80E-12,0.0E+00,-310.0))
  RCONST(247) = (GCARR(5.50E-12,0.0E+00,200.0))
  RCONST(248) = (GCARR(2.40E-12,0.0E+00,40.0))
  RCONST(249) = (GCARR(2.80E-14,0.0E+00,860.0))
  RCONST(250) = (GCARR(8.80E-12,0.0E+00,260.0))
  RCONST(251) = (GCARR(4.90E-11,0.0E+00,0.0))
  RCONST(252) = (GCARR(2.10E-11,0.0E+00,240.0))
  RCONST(253) = (GCARR(1.70E-11,0.0E+00,250.0))
  RCONST(254) = (GCARR(1.60E-11,0.0E+00,0.0))
  RCONST(255) = (GCARR(1.70E-11,0.0E+00,-800.0))
  RCONST(256) = (GCARR(1.80E-11,0.0E+00,-460.0))
  RCONST(257) = (GCARR(1.66E-10,0.0E+00,-7000.0))
  RCONST(258) = (GCARR(2.36E-10,0.0E+00,-6411.0))
  RCONST(259) = (GCARR(8.77E-11,0.0E+00,-4330.0))
  RCONST(260) = (GCJPLPR(4.20E-31,2.4E+00,0.0,2.7E-11,0.0,0.0,0.6,0.0,0.0))
  RCONST(261) = (GCJPLPR(5.40E-31,3.1E+00,0.0,6.5E-12,2.9,0.0,0.6,0.0,0.0))
  RCONST(262) = (GCARR(9.00E-13,0.0E+00,-360.0))
  RCONST(263) = (GCARR(2.00E-12,0.0E+00,-840.0))
  RCONST(264) = (GCARR(1.42E-12,0.0E+00,-1150.0))
  RCONST(265) = (GCARR(2.60E-12,0.0E+00,-1100.0))
  RCONST(266) = (GCARR(1.80E-11,0.0E+00,-600.0))
  RCONST(267) = (GCARR(7.40E-12,0.0E+00,270.0))
  RCONST(268) = (GCARR(6.00E-13,0.0E+00,230.0))
  RCONST(269) = (GCARR(1.40E-12,0.0E+00,600.0))
  RCONST(270) = (GCARR(6.00E-13,0.0E+00,670.0))
  RCONST(271) = (GCARR(1.80E-12,0.0E+00,-250.0))
  RCONST(272) = (GCARR(3.00E-12,0.0E+00,-500.0))
  RCONST(273) = (GCARR(2.40E-12,0.0E+00,-1250.0))
  RCONST(274) = (GCARR(1.20E-12,0.0E+00,-330.0))
  RCONST(275) = (GCARR(1.96E-12,0.0E+00,-1200.0))
  RCONST(276) = (GCARR(2.61E-12,0.0E+00,-944.0))
  RCONST(277) = (GCARR(4.69E-12,0.0E+00,-1134.0))
  RCONST(278) = (GCARR(7.10E-12,0.0E+00,-1270.0))
  RCONST(279) = (GCARR(7.32E-11,0.0E+00,-30.0))
  RCONST(280) = (GCARR(2.30E-11,0.0E+00,-200.0))
  RCONST(281) = (GCARR(1.10E-11,0.0E+00,-980.0))
  RCONST(282) = (GCARR(1.40E-11,0.0E+00,270.0))
  RCONST(283) = (GCARR(3.60E-11,0.0E+00,-375.0))
  RCONST(284) = (GCARR(2.60E-12,0.0E+00,290.0))
  RCONST(285) = (GCARR(6.40E-12,0.0E+00,290.0))
  RCONST(286) = (GCJPLPR(1.80E-31,3.4E+00,0.0,1.50E-11,1.9E0,0.0,0.6,0.0,0.0))
  RCONST(287) = (GCARR(1.00E-12,0.0E+00,-1590.0))
  RCONST(288) = (GCARR(3.00E-11,0.0E+00,-2450.0))
  RCONST(289) = (GCARR(3.50E-13,0.0E+00,-1370.0))
  RCONST(290) = (GCJPLPR(2.20E-33,3.1E+00,0.0,1.8E-10,0.0,0.0,0.6,0.0,0.0))
  RCONST(291) = (GCJPLPR(3.33E-09,3.1E+00,-2502.0,2.73E+14,0.0,-2502.0,0.6,0.0,0.0))
  RCONST(292) = (GCJPLPR(1.90E-32,3.6E+00,0.0,3.7E-12,1.6,0.0,0.6,0.0,0.0))
  RCONST(293) = (GCJPLPR(9.30E-06,4.5E+00,-8649.0,1.74E+15,2.0E0,-8649.0,0.6,0.0,0.0))
  RCONST(294) = (GCARR(2.30E-10,0.0E+00,0.0))
  RCONST(295) = (GCARR(1.20E-11,0.0E+00,0.0))
  RCONST(296) = (GCARR(9.50E-13,0.0E+00,550.0))
  RCONST(297) = (GCARR(2.30E-12,0.0E+00,260.0))
  RCONST(298) = (GCARR(4.10E-13,0.0E+00,290.0))
  RCONST(299) = (GCARR(6.50E-12,0.0E+00,135.0))
  RCONST(300) = (GCARR(2.17E-11,0.0E+00,-1130.0))
  RCONST(301) = (GCARR(1.24E-12,0.0E+00,-1070.0))
  RCONST(302) = (GCARR(3.77E-12,0.0E+00,-1011.0))
  RCONST(303) = (GCARR(2.00E-13,0.0E+00,0.0))
  RCONST(304) = (GCARR(1.60E-10,0.0E+00,0.0))
  RCONST(305) = (GCARR(5.7E-11,0.0E+00,0.0))
  RCONST(306) = (GCARR(7.2E-11,0.0E+00,-70.0))
  RCONST(307) = (GCARR(7.4E-11,0.0E+00,0.0))
  RCONST(308) = (GCARR(5.5E-11,0.0E+00,0.0))
  RCONST(309) = (GCARR(9.6E-11,0.0E+00,0.0))
  RCONST(310) = (GCARR(2.8E-14,0.0E+00,0.0))
  RCONST(311) = (GCARR(6.54E-11,0.0E+00,60.0))
  RCONST(312) = (GCARR(8.12E-11,0.0E+00,-90.0))
  RCONST(313) = (GCARR(7.70E-11,0.0E+00,-1000.0))
  RCONST(314) = (GCARR(7.60E-11,0.0E+00,500.0))
  RCONST(315) = (GCARR(2.05E-10,0.0E+00,0.0))
  RCONST(316) = (GCJPLPR(4.00E-28,0.0E+00,0.0,2.8E-10,0.0E0,0.0,0.6,0.0,0.0))
  RCONST(317) = (GCARR(3.60E-12,0.0E+00,0.0))
  RCONST(318) = (GCJPLPR(1.80E-32,1.0E+00,0.0,1.77E-11,0.0,0.0,0.6,0.0,0.0))
  RCONST(319) = (GCARR(8.40E-11,0.0E+00,-2620.0))
  RCONST(320) = (GCJPLPR(3.00E-31,1.0E+00,0.0,6.6E-11,0.0,0.0,0.63,0.0,0.0))
  RCONST(321) = (GCARR(9.94E+17,0.0E+00,-11859.0))
  RCONST(322) = (GCARR(2.90E-11,0.0E+00,-2600.0))
  RCONST(323) = (GCARR(1.50E-12,0.0E+00,0.0))
  RCONST(324) = (GCJPLPR(7.50E-31,3.5E+00,0.0,7.6E-12,1.5E0,0.0,0.6,0.0,0.0))
  RCONST(325) = (GCARR(2.10E+15,0.0E+00,-13670.0))
  RCONST(326) = (GCARR(9.10E-11,0.0E+00,-146.0))
  RCONST(327) = (GCARR(1.20E-11,0.0E+00,0.0))
  RCONST(328) = (GCARR(3.00E-12,0.0E+00,510.0))
  RCONST(329) = (GCARR(1.20E-11,0.0E+00,510.0))
  RCONST(330) = (GCARR(1.00E-10,0.0E+00,0.0))
  RCONST(331) = (GCARR(1.50E-10,0.0E+00,0.0))
  RCONST(332) = (GCARR(3.80E-02,0.0E+00,0.0))
  RCONST(333) = (GCARR(1.10E-12,0.0E+00,542.0))
  RCONST(334) = (GCARR(5.10E-12,0.0E+00,280.0))
  RCONST(335) = (GCARR(2.81E-12,0.0E+00,280.0))
  RCONST(336) = (GCARR(1.02E-12,0.0E+00,280.0))
  RCONST(337) = (GCARR(2.30E-11,0.0E+00,-870.0))
  RCONST(338) = (GCARR(1.50E-11,0.0E+00,-1090.0))
  RCONST(339) = (GCARR(1.80E-10,0.0E+00,0.0))
  RCONST(340) = (GCARR(3.00E-11,0.0E+00,0.0))
  RCONST(341) = (GCARR(5.00E-12,0.0E+00,0.0))
  RCONST(342) = (GCARR(1.30E-11,0.0E+00,570.0))
  RCONST(343) = (GCARR(9.10E-12,0.0E+00,240.0))
  RCONST(344) = (GCARR(6.00E-12,0.0E+00,500.0))
  RCONST(345) = (GCARR(9.00E-12,0.0E+00,500.0))
  RCONST(346) = (GCARR(1.00E+12,0.0E+00,-9770.0))
  RCONST(347) = (GCARR(2.50E+14,0.0E+00,-9770.0))
  RCONST(348) = (GCARR(2.90E-12,0.0E+00,-1100.0))
  RCONST(349) = (GCARR(3.35E-11,0.0E+00,380.0))
  RCONST(350) = (GC_RO2NO('B',2.70E-12,0.0E+00,350.0,5.0,0.0,0.0))
  RCONST(351) = (GC_RO2NO('A',2.70E-12,0.0E+00,350.0,5.0,0.0,0.0))
  RCONST(352) = (GC_RO2HO2(2.91E-13,0.0E+00,1300.0,5.0,0.0,0.0))
  RCONST(353) = (GCARR(1.20E-11,0.0E+00,652.0))
  RCONST(354) = (GCARR(2.40E-12,0.0E+00,745.0))
  RCONST(355) = (GCARR(2.40E-12,0.0E+00,360.0))
  RCONST(356) = (GCARR(8.70E-14,0.0E+00,1650.0))
  RCONST(357) = (GCARR(2.40E-12,0.0E+00,360.0))
  RCONST(358) = (GCARR(8.70E-14,0.0E+00,1650.0))
  RCONST(359) = (GCARR(4.75E-12,0.0E+00,200.0))
  RCONST(360) = (GCARR(4.40E-13,0.0E+00,380.0))
  RCONST(361) = (GCARR(8.79E-13,0.0E+00,380.0))
  RCONST(362) = (GCARR(2.70E-12,0.0E+00,350.0))
  RCONST(363) = (GC_RO2HO2(2.91E-13,0.0E+00,1300.0,4.0,0.0,0.0))
  RCONST(364) = (GCJPLPR(9.00E-28,8.9E+00,0.0,7.7E-12,0.2,0.0,0.6,0.0,0.0))
  RCONST(365) = (GCJPLEQ(9.00E-29,0.0E+00,14000.0,9.00E-28,8.9E0,0.0,7.7E-12,0.2,0.,0.6,0.,0.))
  RCONST(366) = (GCARR(2.52E-12,0.0E+00,410.0))
  RCONST(367) = (GC_RO2NO('B',2.70E-12,0.0E+00,350.0,5.0,0.0,0.0))
  RCONST(368) = (GC_RO2NO('A',2.70E-12,0.0E+00,350.0,5.0,0.0,0.0))
  RCONST(369) = (GC_RO2HO2(2.91E-13,0.0E+00,1300.0,5.0,0.0,0.0))
  RCONST(370) = (GCARR(2.79E-11,0.0E+00,380.0))
  RCONST(371) = (GC_RO2NO('B',2.70E-12,0.0E+00,350.0,5.0,0.0,0.0))
  RCONST(372) = (GC_RO2NO('A',2.70E-12,0.0E+00,350.0,5.0,0.0,0.0))
  RCONST(373) = (GC_RO2HO2(2.91E-13,0.0E+00,1300.0,4.0,0.0,0.0))
  RCONST(374) = (GCARR(2.00E-17,0.0E+00,0.0))
  RCONST(375) = (GCARR(2.40E-12,0.0E+00,0.0))
  RCONST(376) = (GCARR(6.70E-13,0.0E+00,0.0))
  RCONST(377) = (GCARR(1.20E-15,0.0E+00,0.0))
  RCONST(378) = (GCARR(1.00E-14,0.0E+00,0.0))
  RCONST(379) = (GCARR(1.00E-15,0.0E+00,0.0))
  RCONST(380) = (GCARR(7.00E-14,0.0E+00,0.0))
  RCONST(381) = (GCARR(6.00E-18,0.0E+00,0.0))
  RCONST(382) = (GCARR(1.00E-17,0.0E+00,0.0))
  RCONST(383) = (GCARR(1.20E-15,0.0E+00,0.0))
  RCONST(384) = (GCARR(1.00E-14,0.0E+00,0.0))
  RCONST(385) = (GCARR(1.00E-15,0.0E+00,0.0))
  RCONST(386) = (GCARR(7.00E-14,0.0E+00,0.0))
  RCONST(387) = (GCARR(1.00E-17,0.0E+00,0.0))
  RCONST(388) = (GCARR(6.00E-18,0.0E+00,0.0))
  RCONST(389) = (GCARR(1.20E-15,0.0E+00,0.0))
  RCONST(390) = (GCARR(1.00E-14,0.0E+00,0.0))
  RCONST(391) = (GCARR(1.00E-15,0.0E+00,0.0))
  RCONST(392) = (GCARR(7.00E-14,0.0E+00,0.0))
  RCONST(393) = (GCARR(6.00E-18,0.0E+00,0.0))
  RCONST(394) = (GCARR(1.20E-15,0.0E+00,0.0))
  RCONST(395) = (GCARR(1.00E-14,0.0E+00,0.0))
  RCONST(396) = (GCARR(1.00E-15,0.0E+00,0.0))
  RCONST(397) = (GCARR(7.00E-14,0.0E+00,0.0))
  RCONST(398) = (GCARR(6.00E-18,0.0E+00,0.0))
  RCONST(399) = (GCARR(1.20E-15,0.0E+00,0.0))
  RCONST(400) = (GCARR(1.00E-14,0.0E+00,0.0))
  RCONST(401) = (GCARR(1.00E-15,0.0E+00,0.0))
  RCONST(402) = (GCARR(7.00E-14,0.0E+00,0.0))
  RCONST(403) = (GCARR(1.00E-17,0.0E+00,0.0))
  RCONST(404) = (GCARR(6.00E-18,0.0E+00,0.0))
  RCONST(405) = (GCARR(1.20E-15,0.0E+00,0.0))
  RCONST(406) = (GCARR(1.00E-14,0.0E+00,0.0))
  RCONST(407) = (GCARR(1.00E-15,0.0E+00,0.0))
  RCONST(408) = (GCARR(7.00E-14,0.0E+00,0.0))
  RCONST(409) = (GCARR(6.00E-18,0.0E+00,0.0))
  RCONST(410) = (GCARR(1.00E-17,0.0E+00,0.0))
  RCONST(411) = (GCARR(1.20E-15,0.0E+00,0.0))
  RCONST(412) = (GCARR(1.00E-14,0.0E+00,0.0))
  RCONST(413) = (GCARR(1.00E-15,0.0E+00,0.0))
  RCONST(414) = (GCARR(7.00E-14,0.0E+00,0.0))
  RCONST(415) = (GCARR(6.00E-18,0.0E+00,0.0))
  RCONST(416) = (GCARR(1.00E-17,0.0E+00,0.0))
  RCONST(417) = (GCARR(1.21E-11,0.0E+00,440.0))
  RCONST(418) = (GCARR(1.21E-11,0.0E+00,440.0))
  RCONST(419) = (GCARR(4.00E-12,0.0E+00,0.0))
  RCONST(420) = (GCARR(1.50E-11,0.0E+00,0.0))
  RCONST(421) = (GCARR(3.56E-14,0.0E+00,708.0))
  RCONST(422) = (GCARR(7.40E-13,0.0E+00,765.0))
  RCONST(423) = (GCARR(1.20E-12,0.0E+00,0.0))
  RCONST(424) = (GCARR(5.00E-16,0.0E+00,-530.0))
  RCONST(425) = (GCARR(5.00E-16,0.0E+00,-530.0))
  RCONST(426) = (GCARR(8.33E-13,0.0E+00,490.0))
  RCONST(427) = (GCARR(8.33E-13,0.0E+00,490.0))
  RCONST(428) = (GCARR(4.20E-11,0.0E+00,401.0))
  RCONST(429) = (GCARR(2.95E-15,0.0E+00,-783.0))
  RCONST(430) = (GCARR(1.22E-11,0.0E+00,0.0))
  RCONST(431) = (GCARR(4.00E-12,0.0E+00,0.0))
  RCONST(432) = (GCARR(1.50E-11,0.0E+00,0.0))
  RCONST(433) = (GCARR(3.56E-14,0.0E+00,708.0))
  RCONST(434) = (GCARR(7.40E-13,0.0E+00,765.0))
  RCONST(435) = (GCARR(1.20E-12,0.0E+00,0.0))
  RCONST(436) = (GCARR(3.40E-12,0.0E+00,190.0))
  RCONST(437) = (GCARR(4.00E-12,0.0E+00,0.0))
  RCONST(438) = (GCARR(4.00E-12,0.0E+00,0.0))
  RCONST(439) = (GCARR(1.66E-13,0.0E+00,1300.0))
  RCONST(440) = (GCARR(1.66E-13,0.0E+00,1300.0))
  RCONST(441) = (GCARR(1.60E-13,0.0E+00,708.0))
  RCONST(442) = (GCARR(9.68E-14,0.0E+00,708.0))
  RCONST(443) = (GCARR(8.85E-13,0.0E+00,765.0))
  RCONST(444) = (GCARR(5.37E-13,0.0E+00,765.0))
  RCONST(445) = (GCARR(1.20E-12,0.0E+00,0.0))
  RCONST(446) = (GCARR(1.20E-12,0.0E+00,0.0))
  RCONST(447) = (GCARR(7.00E-14,0.0E+00,1000.0))
  RCONST(448) = (GCARR(4.25E-14,0.0E+00,1000.0))
  RCONST(449) = (GCARR(2.96E-14,0.0E+00,1000.0))
  RCONST(450) = (GCARR(4.80E-12,0.0E+00,0.0))
  RCONST(451) = (GCARR(7.29E-11,0.0E+00,0.0))
  RCONST(452) = (GCARR(1.67E-16,0.0E+00,0.0))
  RCONST(453) = (GCARR(3.15E-13,0.0E+00,-448.0))
  RCONST(454) = (GCARR(3.15E-13,0.0E+00,-448.0))
  RCONST(455) = (GCARR(2.78E-04,0.0E+00,0.0))
  RCONST(456) = (GCARR(2.78E-04,0.0E+00,0.0))
  RCONST(457) = (GC_OHHNO3(2.41E-14,0.0E+00,460.0,2.69E-17,0.E0,2199.,6.51E-34,0.E0,1335.0))
  RCONST(458) = (GCARR(1.40E-11,0.0E+00,0.0))
  RCONST(459) = (GCARR(5.10E-11,0.0E+00,0.0))
  RCONST(460) = (GCARR(2.70E-12,0.0E+00,360.0))
  RCONST(461) = (GCARR(2.70E-12,0.0E+00,360.0))
  RCONST(462) = (GCARR(2.05E-13,0.0E+00,1300.0))
  RCONST(463) = (GCARR(2.05E-13,0.0E+00,1300.0))
  RCONST(464) = (GCARR(2.90E+07,0.0E+00,-5300.0))
  RCONST(465) = (GCARR(1.28E+07,0.0E+00,-5300.0))
  RCONST(466) = (HET(ind_HO2,1))
  RCONST(467) = (HET(ind_NO2,1))
  RCONST(468) = (HET(ind_NO3,1))
  RCONST(469) = (HET(ind_N2O5,1))
  RCONST(470) = (HET(ind_N2O5,2))
  RCONST(471) = (HET(ind_N2O5,3))
  RCONST(472) = (HET(ind_BrNO3,1))
  RCONST(473) = (HET(ind_BrNO3,2))
  RCONST(474) = (HET(ind_ClNO3,1))
  RCONST(475) = (HET(ind_ClNO3,2))
  RCONST(476) = (HET(ind_ClNO3,3))
  RCONST(477) = (HET(ind_ClNO3,4))
  RCONST(478) = (HET(ind_ClNO3,5))
  RCONST(479) = (HET(ind_HOCl,1))
  RCONST(480) = (HET(ind_HOCl,2))
  RCONST(481) = (HET(ind_HOBr,1))
  RCONST(482) = (HET(ind_HOBr,2))
  RCONST(483) = (HET(ind_HOBr,3))
  RCONST(484) = (HET(ind_HOBr,4))
  RCONST(485) = (HET(ind_HOBr,5))
  RCONST(486) = (HET(ind_HOBr,6))
  RCONST(487) = (HET(ind_HOBr,7))
  RCONST(488) = (HET(ind_HCl,1))
  RCONST(489) = (HET(ind_HCl,2))
  RCONST(490) = (HET(ind_O3,1))
  RCONST(491) = (HET(ind_O3,2))
  RCONST(492) = (HET(ind_O3,3))
  RCONST(493) = (HET(ind_HBr,1))
  RCONST(494) = (HET(ind_HBr,2))
  RCONST(495) = (HET(ind_HI,1))
  RCONST(496) = (HET(ind_HI,2))
  RCONST(497) = (HET(ind_HI,3))
  RCONST(498) = (HET(ind_HOI,1))
  RCONST(499) = (HET(ind_HOI,2))
  RCONST(500) = (HET(ind_I2O2,1))
  RCONST(501) = (HET(ind_I2O2,2))
  RCONST(502) = (HET(ind_I2O2,3))
  RCONST(503) = (HET(ind_I2O3,1))
  RCONST(504) = (HET(ind_I2O3,2))
  RCONST(505) = (HET(ind_I2O3,3))
  RCONST(506) = (HET(ind_I2O4,1))
  RCONST(507) = (HET(ind_I2O4,2))
  RCONST(508) = (HET(ind_I2O4,3))
  RCONST(509) = (HET(ind_IONO,1))
  RCONST(510) = (HET(ind_IONO,2))
  RCONST(511) = (HET(ind_IONO2,1))
  RCONST(512) = (HET(ind_IONO2,2))
  RCONST(513) = (HET(ind_IONO,3))
  RCONST(514) = (HET(ind_IONO2,3))
  RCONST(515) = (HET(ind_HOI,3))
  RCONST(516) = (HET(ind_GLYX,1))
  RCONST(517) = (HET(ind_MGLY,1))
  RCONST(518) = (HET(ind_IEPOXA,1))
  RCONST(519) = (HET(ind_IEPOXB,1))
  RCONST(520) = (HET(ind_IEPOXD,1))
  RCONST(521) = (HET(ind_IMAE,1))
  RCONST(522) = (HET(ind_LVOC,1))
  RCONST(523) = (HET(ind_ISN1OG,1))
  RCONST(524) = (HET(ind_ISOPND,1))
  RCONST(525) = (HET(ind_ISOPNB,1))
  RCONST(526) = (HET(ind_MACRN,1))
  RCONST(527) = (HET(ind_MVKN,1))
  RCONST(528) = (HET(ind_R4N2,1))
  RCONST(529) = (HET(ind_DHDN,1))
  RCONST(530) = (HET(ind_MONITS,1))
  RCONST(531) = (HET(ind_MONITU,1))
  RCONST(532) = (HET(ind_HONIT,1))
  RCONST(533) = (PHOTOL(3))
  RCONST(534) = (PHOTOL(103))
  RCONST(535) = (PHOTOL(11))
  RCONST(536) = (PHOTOL(9))
  RCONST(537) = (PHOTOL(10))
  RCONST(538) = (PHOTOL(7))
  RCONST(539) = (PHOTOL(8))
  RCONST(540) = (PHOTOL(16))
  RCONST(541) = (PHOTOL(15))
  RCONST(542) = (PHOTOL(17))
  RCONST(543) = (PHOTOL(18))
  RCONST(544) = (PHOTOL(12))
  RCONST(545) = (PHOTOL(13))
  RCONST(546) = (PHOTOL(14))
  RCONST(547) = (PHOTOL(61))
  RCONST(548) = (PHOTOL(62))
  RCONST(549) = (PHOTOL(59))
  RCONST(550) = (PHOTOL(70))
  RCONST(551) = (PHOTOL(76))
  RCONST(552) = (PHOTOL(77))
  RCONST(553) = (PHOTOL(69))
  RCONST(554) = (PHOTOL(68))
  RCONST(555) = (PHOTOL(72))
  RCONST(556) = (PHOTOL(73))
  RCONST(557) = (PHOTOL(74))
  RCONST(558) = (PHOTOL(71))
  RCONST(559) = (PHOTOL(63))
  RCONST(560) = (PHOTOL(64))
  RCONST(561) = (PHOTOL(65))
  RCONST(562) = (PHOTOL(66))
  RCONST(563) = (PHOTOL(75))
  RCONST(564) = (PHOTOL(78))
  RCONST(565) = (PHOTOL(79))
  RCONST(566) = (PHOTOL(80))
  RCONST(567) = (PHOTOL(81))
  RCONST(568) = (PHOTOL(82))
  RCONST(569) = (PHOTOL(83))
  RCONST(570) = (PHOTOL(84))
  RCONST(571) = (PHOTOL(85))
  RCONST(572) = (PHOTOL(86))
  RCONST(573) = (PHOTOL(86))
  RCONST(574) = (PHOTOL(86))
  RCONST(575) = (PHOTOL(87))
  RCONST(576) = (PHOTOL(88))
  RCONST(577) = (PHOTOL(89))
  RCONST(578) = (PHOTOL(90))
  RCONST(579) = (PHOTOL(91))
  RCONST(580) = (PHOTOL(98))
  RCONST(581) = (PHOTOL(99))
  RCONST(582) = (PHOTOL(92))
  RCONST(583) = (PHOTOL(93))
  RCONST(584) = (PHOTOL(94))
  RCONST(585) = (PHOTOL(23))
  RCONST(586) = (PHOTOL(28))
  RCONST(587) = (PHOTOL(32))
  RCONST(588) = (PHOTOL(29))
  RCONST(589) = (PHOTOL(30))
  RCONST(590) = (PHOTOL(31))
  RCONST(591) = (PHOTOL(56))
  RCONST(592) = (PHOTOL(55))
  RCONST(593) = (PHOTOL(50))
  RCONST(594) = (PHOTOL(43))
  RCONST(595) = (PHOTOL(45))
  RCONST(596) = (PHOTOL(33))
  RCONST(597) = (PHOTOL(22))
  RCONST(598) = (PHOTOL(27))
  RCONST(599) = (PHOTOL(25))
  RCONST(600) = (PHOTOL(26))
  RCONST(601) = (PHOTOL(21))
  RCONST(602) = (PHOTOL(19))
  RCONST(603) = (PHOTOL(20))
  RCONST(604) = (PHOTOL(24))
  RCONST(605) = (PHOTOL(102))
  RCONST(606) = (PHOTOL(115))
  RCONST(607) = (PHOTOL(116))
  RCONST(608) = (PHOTOL(117))
  RCONST(609) = (PHOTOL(118))
  RCONST(610) = (PHOTOL(119))
  RCONST(611) = (PHOTOL(120))
  RCONST(612) = (PHOTOL(121))
  RCONST(613) = (PHOTOL(122))
  RCONST(614) = (PHOTOL(123))
  RCONST(615) = (PHOTOL(124))
  RCONST(616) = (PHOTOL(125))
  RCONST(617) = (PHOTOL(126))
  RCONST(618) = (PHOTOL(127))
  RCONST(619) = (PHOTOL(128))
  RCONST(620) = (PHOTOL(129))
  RCONST(621) = (PHOTOL(130))
  RCONST(622) = (PHOTOL(104))
  RCONST(623) = (PHOTOL(105))
  RCONST(624) = (PHOTOL(95))
  RCONST(625) = (PHOTOL(96))
  RCONST(626) = (PHOTOL(97))
  RCONST(627) = (PHOTOL(106))
  RCONST(628) = (PHOTOL(107))
  RCONST(629) = (PHOTOL(108))
  RCONST(630) = (PHOTOL(109))
  RCONST(631) = (PHOTOL(110))
  RCONST(632) = (PHOTOL(111))
  RCONST(633) = (PHOTOL(112))
  RCONST(634) = (PHOTOL(113))
  RCONST(635) = (PHOTOL(114))
  RCONST(636) = (GCARR(9.10E-15,0.0E+00,-2580.0))
  RCONST(637) = (GCJPLPR(1.10E-28,3.5E+00,0.0,8.4E-12,1.75,0.0,0.5,0.0,0.0))
  RCONST(638) = (GCARR(1.53E-13,0.0E+00,1300.0))
  RCONST(639) = (GCARR(2.7E-12,0.0E+00,360.0))
! RCONST(640) = constant rate coefficient
! RCONST(641) = constant rate coefficient
  RCONST(642) = (GCARR(9.5E-13,0.0E+00,-5988.0))
  RCONST(643) = (GCARR(2.5E-14,0.0E+00,-300.0))
! RCONST(644) = constant rate coefficient
  RCONST(645) = (GCARR(1.90E-12,0.0E+00,190.0))
! RCONST(646) = constant rate coefficient
  RCONST(647) = (GCJPLPR(5.50E-30,0.0E+00,0.0,8.3E-13,-2.0,0.0,0.5,0.0,0.0))
  RCONST(648) = (PHOTOL(134))
  RCONST(649) = (GCARR(2.3E-12,0.0E+00,-193.0))
  RCONST(650) = (GCARR(1.8E-12,0.0E+00,340.0))
! RCONST(651) = constant rate coefficient
  RCONST(652) = (2.91D-13*EXP(1300/TEMP)*0.82)
  RCONST(653) = (GCARR(2.60E-12,0.0E+00,365.0))
! RCONST(654) = constant rate coefficient
  RCONST(655) = (GCARR(1.70E-14,0.0E+00,220.0))
  RCONST(656) = (GCARR(4.20E-14,0.0E+00,220.0))
  RCONST(657) = (GCARR(4.70E-13,0.0E+00,1220.0))
! RCONST(658) = constant rate coefficient
! RCONST(659) = constant rate coefficient
! RCONST(660) = constant rate coefficient
! RCONST(661) = constant rate coefficient
! RCONST(662) = constant rate coefficient
! RCONST(663) = constant rate coefficient
  RCONST(664) = (GCARR(5.90E-12,0.0E+00,225.0))
! RCONST(665) = constant rate coefficient
  RCONST(666) = (GCARR(1.10E-11,0.0E+00,340.0))
  RCONST(667) = (GCARR(7.50E-12,0.0E+00,290.0))
  RCONST(668) = (PAN(3.28E-28,0.,-6.87,1.125E-11,0.,-1.105,0.3))
! RCONST(669) = constant rate coefficient
  RCONST(670) = (PAN(1.10E-5,-10100.,0.,1.90E+17,-14100.,0.,0.3)*0.67)
! RCONST(671) = constant rate coefficient
! RCONST(672) = constant rate coefficient
  RCONST(673) = (GCARR(2.670E-12,0.0E+00,365.0))
! RCONST(674) = constant rate coefficient
  RCONST(675) = (GCARR(2.24E-13,0.0E+00,1300.0))
! RCONST(676) = constant rate coefficient
! RCONST(677) = constant rate coefficient
! RCONST(678) = constant rate coefficient
! RCONST(679) = constant rate coefficient
! RCONST(680) = constant rate coefficient
  RCONST(681) = (GCARR(2.670E-13,0.0E+00,365.0))
  RCONST(682) = (GCARR(2.670E-12,0.0E+00,365.0))
! RCONST(683) = constant rate coefficient
! RCONST(684) = constant rate coefficient
! RCONST(685) = constant rate coefficient
! RCONST(686) = constant rate coefficient
! RCONST(687) = constant rate coefficient
! RCONST(688) = constant rate coefficient
  RCONST(689) = (PHOTOL(133))
  RCONST(690) = (PHOTOL(132))
  RCONST(691) = (PHOTOL(131))
  RCONST(692) = (PHOTOL(135))
  RCONST(693) = (HET(ind_MONITU,1))
! RCONST(694) = constant rate coefficient
! RCONST(695) = constant rate coefficient
! RCONST(696) = constant rate coefficient
! RCONST(697) = constant rate coefficient
! RCONST(698) = constant rate coefficient
! RCONST(699) = constant rate coefficient
      
END SUBROUTINE Update_RCONST

! End of Update_RCONST function
! ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~


! ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
! 
! Update_PHOTO - function to update photolytical rate constants
!   Arguments :
! 
! ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~

SUBROUTINE Update_PHOTO ( )


   USE gckpp_Global

  RCONST(533) = (PHOTOL(3))
  RCONST(534) = (PHOTOL(103))
  RCONST(535) = (PHOTOL(11))
  RCONST(536) = (PHOTOL(9))
  RCONST(537) = (PHOTOL(10))
  RCONST(538) = (PHOTOL(7))
  RCONST(539) = (PHOTOL(8))
  RCONST(540) = (PHOTOL(16))
  RCONST(541) = (PHOTOL(15))
  RCONST(542) = (PHOTOL(17))
  RCONST(543) = (PHOTOL(18))
  RCONST(544) = (PHOTOL(12))
  RCONST(545) = (PHOTOL(13))
  RCONST(546) = (PHOTOL(14))
  RCONST(547) = (PHOTOL(61))
  RCONST(548) = (PHOTOL(62))
  RCONST(549) = (PHOTOL(59))
  RCONST(550) = (PHOTOL(70))
  RCONST(551) = (PHOTOL(76))
  RCONST(552) = (PHOTOL(77))
  RCONST(553) = (PHOTOL(69))
  RCONST(554) = (PHOTOL(68))
  RCONST(555) = (PHOTOL(72))
  RCONST(556) = (PHOTOL(73))
  RCONST(557) = (PHOTOL(74))
  RCONST(558) = (PHOTOL(71))
  RCONST(559) = (PHOTOL(63))
  RCONST(560) = (PHOTOL(64))
  RCONST(561) = (PHOTOL(65))
  RCONST(562) = (PHOTOL(66))
  RCONST(563) = (PHOTOL(75))
  RCONST(564) = (PHOTOL(78))
  RCONST(565) = (PHOTOL(79))
  RCONST(566) = (PHOTOL(80))
  RCONST(567) = (PHOTOL(81))
  RCONST(568) = (PHOTOL(82))
  RCONST(569) = (PHOTOL(83))
  RCONST(570) = (PHOTOL(84))
  RCONST(571) = (PHOTOL(85))
  RCONST(572) = (PHOTOL(86))
  RCONST(573) = (PHOTOL(86))
  RCONST(574) = (PHOTOL(86))
  RCONST(575) = (PHOTOL(87))
  RCONST(576) = (PHOTOL(88))
  RCONST(577) = (PHOTOL(89))
  RCONST(578) = (PHOTOL(90))
  RCONST(579) = (PHOTOL(91))
  RCONST(580) = (PHOTOL(98))
  RCONST(581) = (PHOTOL(99))
  RCONST(582) = (PHOTOL(92))
  RCONST(583) = (PHOTOL(93))
  RCONST(584) = (PHOTOL(94))
  RCONST(585) = (PHOTOL(23))
  RCONST(586) = (PHOTOL(28))
  RCONST(587) = (PHOTOL(32))
  RCONST(588) = (PHOTOL(29))
  RCONST(589) = (PHOTOL(30))
  RCONST(590) = (PHOTOL(31))
  RCONST(591) = (PHOTOL(56))
  RCONST(592) = (PHOTOL(55))
  RCONST(593) = (PHOTOL(50))
  RCONST(594) = (PHOTOL(43))
  RCONST(595) = (PHOTOL(45))
  RCONST(596) = (PHOTOL(33))
  RCONST(597) = (PHOTOL(22))
  RCONST(598) = (PHOTOL(27))
  RCONST(599) = (PHOTOL(25))
  RCONST(600) = (PHOTOL(26))
  RCONST(601) = (PHOTOL(21))
  RCONST(602) = (PHOTOL(19))
  RCONST(603) = (PHOTOL(20))
  RCONST(604) = (PHOTOL(24))
  RCONST(605) = (PHOTOL(102))
  RCONST(606) = (PHOTOL(115))
  RCONST(607) = (PHOTOL(116))
  RCONST(608) = (PHOTOL(117))
  RCONST(609) = (PHOTOL(118))
  RCONST(610) = (PHOTOL(119))
  RCONST(611) = (PHOTOL(120))
  RCONST(612) = (PHOTOL(121))
  RCONST(613) = (PHOTOL(122))
  RCONST(614) = (PHOTOL(123))
  RCONST(615) = (PHOTOL(124))
  RCONST(616) = (PHOTOL(125))
  RCONST(617) = (PHOTOL(126))
  RCONST(618) = (PHOTOL(127))
  RCONST(619) = (PHOTOL(128))
  RCONST(620) = (PHOTOL(129))
  RCONST(621) = (PHOTOL(130))
  RCONST(622) = (PHOTOL(104))
  RCONST(623) = (PHOTOL(105))
  RCONST(624) = (PHOTOL(95))
  RCONST(625) = (PHOTOL(96))
  RCONST(626) = (PHOTOL(97))
  RCONST(627) = (PHOTOL(106))
  RCONST(628) = (PHOTOL(107))
  RCONST(629) = (PHOTOL(108))
  RCONST(630) = (PHOTOL(109))
  RCONST(631) = (PHOTOL(110))
  RCONST(632) = (PHOTOL(111))
  RCONST(633) = (PHOTOL(112))
  RCONST(634) = (PHOTOL(113))
  RCONST(635) = (PHOTOL(114))
  RCONST(648) = (PHOTOL(134))
  RCONST(689) = (PHOTOL(133))
  RCONST(690) = (PHOTOL(132))
  RCONST(691) = (PHOTOL(131))
  RCONST(692) = (PHOTOL(135))
      
END SUBROUTINE Update_PHOTO

! End of Update_PHOTO function
! ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~



END MODULE gckpp_Rates

